The Lagrangian description of perfect fluids and modified gravity with an extra force 
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We revisit the issue of the correct Lagrangian description of a perfect fluid (£i = P versus 
£2 = —p) in relation with modified gravity theories in which galactic luminous matter couples 
nonminimally to the Ricci scalar. These Lagrangians are only equivalent when the fluid couples 
minimally to gravity and not otherwise; in the presence of nonminimal coupling they give rise to 
two distinct theories with different predictions. 
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INTRODUCTION 

Attempts to explain the present acceleration of the uni- 
verse discovered with type la supernovae [l| without in- 
voking a mysterious and exotic dark energy Q have lead 
to modifying gravity at the largest scales in the so-called 
f(R) gravity (0], see [i] for a comprehensive review and 
[5| for brief introductions). f{R) theories of gravity have 
been employed also in attempts to explain dark matter in 
galaxies and clusters Q. A recent proposal ny Bertolami, 
Bdhmer, Harko, and Lobo (BBHL) [Tj contemplates the 
possibility of coupling matter explicitly to the Ricci cur- 
vature, as described by the action [l^ 
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where R is the Ricci curvature of spacetime, fi^2{R) are 
functions of i?, A is a small parameter, and is the 
matter Lagrangian density (see Q for similar theories 
and [13] for the special case of a gauge field coupling 
explicitly to the curvature). The extra coupling leads to 
the non-conservation of the matter stress-energy tensor 
^ir^ , according to 
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where a prime denotes differentiation with respect to R. 
In [3] BBHL found that an extra force, not appearing in 
the minimally coupled versions of f{R) gravity, acts upon 
particles and could effectively replace dark matter. This 
intriguing possibility has motivated further work ll|, |12| . 

n. 

In 0] a dust fluid model is assumed to describe ordi- 
nary (luminous) matter in galaxies and the extra force 



P = adequately reproducing the non-relativistic mo- 
tions of stars, making again dark matter a necessary in- 
gredient to explain the galactic rotation curves. Sub- 
sequently Bertolami and collaborators pointed out that 
an equivalent Lagrangian density for a perfect fluid is 
£2 = —p, which is obtained by adding surface terms to 
the action = / d^x \J—gL\ [li!] (see for detailed 
studies of the Lagrangian formalism for perfect fluids 
(20|). From [12] it would appear that the extra force 
that has the potential to replace dark matter is present 
or absent according to the choice that is made between 
the two supposedly equivalent Lagrangians. This shows 
that, clearly, the two Lagrangian densities Ci = P and 
C2 — —p cannot be equivalent, and the recent literature 
has discussed this issue in terms of the problem of "which 
Lagrangian density correctly describes a perfect fluid" . 
Here we show that, posed in these terms, this problem is 
meaningless. In fact, for a perfect fluid that docs not cou- 
ple explicitly to the curvature (i.e., for A = 0), the two 
Lagrangian densities Ci — P and £2 — —p a re p erfectly 
equivalent, as discussed extensively in Refs. [l5| and re- 
marked in [13]. However, there is little doubt that for a 
coupled perfect fluid (A ^ 0) the two Lagrangians are in- 
equivalent. This fact is a manifestation of a more general 
situation: if a Lagrangian system consists of two subsys- 
tems and there are two Lagrangians for one of these sub- 
systems, which provide equivalent descriptions for that 
subsystems when it is isolated, the two Lagrangians cease 
to be equivalent when the subsystem couples to the rest 
of the system. An elementary example is provided in the 
next section. We follow the notations of Ref . [id] . 
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TWO COUPLED OSCILLATORS 



was derived (here u'^ is the fluid four-velocity and hab = 
gab + UaUb). In [3] it was noted that, adopting the stan- 
dard Lagrangian density Ci — P for a perfect fluid, this 
extra force vanishes for dust, which has equation of state 



Point particle mechanics provides a very simple ex- 
ample of this situation. Considers two coupled one- 
dimensional oscillators which, in isolation, are described 
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by the Lagrangians 
Li (<7i,gi) 

L2 (92,172) 
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where gi^2 are the respective Lagrangian coordinates 
Now couple the two oscillators via the term L12 = 
—XqiL2 in the total Lagrangian 

1 r 



L (91,^1,92,92) =^1 + ^2 + -^12 = 
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The Euler-Lagrange equations 



yield the equations of motion 
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= 0, 
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92 (1 - A91) - A9192 + K2{1~- A91) 92 = . (9) 

Now consider the new Lagrangian L2 + C (where C is a 
nonzero constant) for the subsystem 2: when this subsys- 
tem is isolated, this new Lagrangian is trivially equivalent 
to L2. However, see what happens when L2 is replaced 
by L2 + C in the description of the coupled subsystem; 
the equations of motion become then 



91 + Kiqi 



(92) - -ft^292 



AC = , 



(10) 



92 (1 - A91) - A9192 + ^^2 (1 - A91) 92 = . (11) 

Eq. pO|) does not coincide with ([8|) . 

To give another, less trivial, example assume that L2 is 
replaced by the equivalent Lagrangian (when the subsys- 
tem 2 is isolated) L2+dF/dt, where F{q2,t) is a function 
of the coordinate 92 and the time t. Then, the total La- 
grangian of the two coupled subsytems becomes 



i* =Li+i2 



(12) 



--A91L2-A91-: 

the change, which leaves the equation of motion for the 
oscillator 2 unaffected when the latter is isolated, now 
affects the way this oscillator interacts with the first one. 
The new equations of motion are 



91 + Ai9i + - 
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+ iv:2(l- A9i)92 =0. 
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Another argument advocated in IJ] as evidence for the 
extra force in galaxies needs to be addressed: in Ref. [13], 
Puetzfeld and Obukhov derive the equations of motion 
of a test body described by a Lagrangian density C and 
an energy-momentum tensor T^™"* in the modified grav- 
ity theories introduced in ,7], using a multipole method. 
They show that, in general, extra forces occur even on 
single-pole particles, as described by eq. (32) of 

[mu' (1 + A/2)] = iV^'AVb/a . (15) 

Here r is the proper time along the worldlinc of a particle 
with four-tangent m° and mass m, while 
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(16) 



S(t) 



where a tilde denotes the density of the corresponding 
quantity and !B°^ = £g°'^ [17[. In spite of being a sig- 



nificant piece of work describing the motion of extended 
objects, the analysis of [l3| does not solve the issue of 
which Lagrangian (£1 or £2) is appropriate to describe 
a perfect fluid: it merely says what are the equations 
of motion once a Lagrangian is assumed. If we assume 
the Lagrangian Ci = P, then S"^ = Pg"* and, conse- 
quently, N"''', vanish for a dust with P — 0. If instead 
the Lagrangian £2 = —p is assumed, N"-'' does not vanish 
and there are extra forces even on the particles of a dust 
fluid. Hence, the work of [l7| des not resolve the issue 
of whether Ci or £2 is appropriate to describe a perfect 
fluid nonminimally coupled to gravity. 

We note that the case of the cosmological constant, re- 
garded as an effective form of matter, is very special: if 
one considers a cosmological constant as a perfect fluid 
with Lagrangian C = —A, there is no difference between 
the two Lagrangians Ci = P and L2 = —p because of the 
peculiar equation of state P — ~p of the effective fluid, 
and there is no extra force on the cosmological constant 
"fluid". This coincidence is, of course, consistent with 
the fact that one can also regard the cosmological con- 
stant term as a purely geometrical term which is part of 
the first function fi{R) in the action ([T]) (for example, 
h{R)^R-K). 

To summarize our conclusions, there is little doubt that 
the two Lagrangian densities Ci = P and £2 = —p are 
equivalent for the description of a perfect fluid which is 
not coupled directly to gravity, as shown in Refs. [l5| . 
However, as soon as this fluid is coupled explicitly to 
gravity as in eq. ([T]), the two Lagrangian densities cease to 
be equivalent. It is not clear that one should be physically 
preferred with respect to the other: simply, they give 
rise to two inequivalent theories of gravity and matter, 
which are both correct. Which one should be chosen 
must be decided by independent arguments. It is a fact 
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that by choosing C = P there is no extra force on a 
dust fluid and it is equally undeniable that by choosing 
£2 = there will be such a force, which may ultimately 
provide an alternative to dark matter. We are not able 
to provide independent arguments in favour of one choice 
or the other: the contribution of the present note merely 
consists in showing that, in spite of the appearance, there 
is no contradiction between the results of UM and those 



of [12| 
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